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Abstract 



"The lateral distribution function (LDF) for coherent electromagnetic radiation from air showers initiated by ultra-high- 
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energy cosmic rays is calculated using a macroscopic description. A new expression is derived to calculate the coherent 
f^'radio pulse at small distances from the observer. It is shown that for small distances to the shower axis the shape of the 
^ electric pulse is determined by the 'pancake' function, describing the longitudinal distribution of charged particles within 
|the shower front, while for large distances the pulse is determined by the shower profile. This reflects in a different scaling 
^\ of the LDF at small and at large distances. As a first application we calculate the LDF for proton- and iron-induced 
T-H showers and we show that this offers a very sensitive measure to discriminate between these two. We show that due to 
■interference between the geo-magnetic and the charge-excess contributions the intensity pattern of the radiation is not 
circular symmetric. 

.Key words: Radio detection. Air showers. Cosmic rays, Geo-synchrotron, Geo-magnetic, Coherent radio emission, 
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i. Introduction 

In recent years much progress has been made modeling 
electric pulses initiated by extensive air showers (EAS) 
13, One of the reasons for this progress are the re- 
sults from the LOPES [1 Q and CODALEMA 0, 0] ex- 
Ipcrimcnts. Both experiments indicate that the dominant 
emission mechanism is due to induction effects from the 
^Earth's magnetic field which exerts a Lorcntz force on 
■the charged particles in the shower. The emission pro- 
cess can be described in a microscopic model where the 
individual electrons and positrons move on cyclotron or- 
bits HBSli^. The importance of coherent radio emis- 
sion was already noted in earlier research on radio emis- 
sion from air showers [ll|, 12, 13, Recently, this ap- 



proach has received renewed attention with the realistic 
calculations in the Macroscopic Geo-Magnctic Radiation 
(MGMR) model as presented in Rcf. . 

The experimental results have triggered plans for an 
extensive array of radio detectors at the Pierre Auger Ob- 
servatory [l5, 3 121- A clear theoretical understanding of 
the pulse shape and its dependence on the distance from 
the shower axis is therefore of importance. In the present 
paper we derive a new expression for the electromagnetic 
pulse that is more appropriate for observer positions that 
are closer to the shower axis while the expression derived 
in Ref. Q is appropriate for large distances. This allows us 
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to address quantitatively the differences in the structure 
of the electromagnetic pulse emitted by the EAS at large 
and small distances. It was already noticed in an earlier 
study 01 that in a macroscopic calculation the structure 
of the electromagnetic pulse at small distances is strongly 
affected by the distribution of the particles in the shower 
front (the pancake). In Ref. [l^] the importance of length 
scales was emphasized in the understanding of macroscopic 
geo-magnetic radiation from the EAS indicating that the 
physics at short distances differs from that at large dis- 
tances. 

In Ref. [l[ it is argued that in the synchrotron-emission 
model REAS2, the lateral distribution function (LDF) can 
be used to disentangle the chemical composition of cos- 
mic rays. As shown in Ref. [3] the predictions of the 
macroscopic model and the synchrotron-emission models 
such as used in Refs. 0, S Hi differ greatly. In par- 
ticular the latter models predict a unipolar pulse which 
contradicts the fact that the intensity of the emitted radi- 
ation should vanish at the longest wavelength a con- 
dition that is satisfied by the bi-polar pulse in the MGMR 
model. Recently 0, UM^ has been shown that this is 
due to the omission of the bremsstrahlung contributions 
at the beginning and the end of the particle trajectories 
in the synchrotron-emission models. This effect is now in- 
cluded in an updated version of the geo-synchrotron code 
REAS3 [i^. Therefore, we will investigate here, as an 
application of the new calculation scheme for the MGMR 
model, the differences in the LDF for iron- and proton- 



Preprint submitted to Elsevier 



August 20, 2010 



induced showers. 

In Section [5] a quick overview of the MGMR model is 
given. For completeness we will shortly review the deriva- 
tion of the expression for the electric field. As the next step 
we concentrate on the short distance scales in Section [3] 
where a new expression is derived for the electric field at 
small impact parameters. For simplicity, we will limit our- 
selves to vertical incoming air showers. The geometry con- 
sidered is an incoming shower with velocity (3 = —jSz, an 
observer placed at a distance d = dx from the shower axis 
with a magnetic field B ~ By perpendicular to the air 
shower. The realistic case will be treated in a forthcoming 
publication. As a first application of the new calculation 
scheme it is shown in Section 2] that the LDF depends on 
the orientation of the observer with respect to the shower 
axis. This angular dependence is due to interference of 
the leading magnetic contribution with secondary contri- 
butions. Of these secondary contributions, the one which 
is generated by charge excess in the shower is the most 
important. As a second application the influence of the 
chemical composition of the cosmic ray on the LDF is dis- 
cussed in Sections 5 and 6. 

2. The model 

When an ultra-high energy cosmic ray collides in the 
atmosphere, a cascade of secondary particles is created, 
moving towards the Earth with a velocity close to that of 
the speed of light. This can be visualized by a 'pancake' of 
particles. The basic picture in the MGMR model [H is that 
the charged particles in the pancake, mostly electrons and 
positrons, will be deflected in the Earth's magnetic fleld 
causing the flow of a macroscopic electric current. The 
strength of this electric current is time dependent and thus 
radiates. We will temporarily make the assumption that 
there is an equal amount of electrons and positrons in the 
shower; in a later section the effects of charge excess will 
be discussed. 

The total amount of electrons and positrons traveling 
in the shower front can be described as a function of height 
z = —cPst' + h, where t' is the shower time (negative) and 
the front of the shower hits the Earth at t' ~ 0. This 
implies that the front of the shower is located at a height of 
z — ~cf3st' with the particles lagging behind at a distance 
within the pancake following a distribution fp{h). The 
distribution of particles is given as 



N{z,t')=NJt{t')fpih) 



(1) 



where ft{t') is the normalized shower profile, and the to- 
tal number of particles at the shower maximum equals 
A^e = 6 X {Ep /lO^^eY) defined b y th e maximum number 
of particles for a 10^^ eV shower j24| . 



The longitudinal profile N^ftit') is parametrized [25| as 
a function of the penetration depth X in units of g cm~^. 



x„ 



-1.5Xlns)/Xo 



(2) 



The penetration depth is written as a function of height 

Using p{0) « 1168 g m-3, and 

-4 ^-1 



as X{z) = {p{0)/C)e 

X{0) « 1000 g cm-2 gives C = 1.168- lO""^ m^^. The pa- 
rameter Xmax is taken to reproduce the shower maximum 
from simulations H], X^ax = (840-^70 logio(i^p/10^° eV)) 
g cm~^. The shower age s{X) can be parameterized as [25| , 

g^^-j ^ 3A/Ao 

A/ Aq + 2Xmax/ Aq 

where Aq = 36.7 g cm~^ is the radiation length of electrons 
in air. 

For the pancake thickness a parametrization can be 
made using measured arrival time distributions. This can 
be fitted using a F-probability distribution function (F- 
pdf) HI, [2^, 

fp{h) = hf'e-^^/'^ X (4/L2), (4) 

where /? = 1 and L ~ 3.9m have been used, see Section[5l 
The electrons and positrons will be deflected in the 
Earth's magnetic field which macroscopically induces a net 
current in the x direction given by, 

fix, y, z, t') = (vaq) e N{z, t')5{x)5{v) , (5) 

where q is the sign of the electric charge e. We disregard 
the lateral distribution of the charged particles and fix the 
position of the charged particles at the shower axis. The 
drift velocity, {vaq) , depends rather strongly on the model 
assumptions made Q and we adopt a value of {vaq) = 
0.025 c. 

The vector potential is given by the Lienard-Wiechert 
fields [i^l in terms of the current density, f , 

ft{tr)f{h) 



A^{t,d)^Jo 



V 



■ dh , 



(6) 



where Jq — {vdq)Nee/ [Aneoc], d is the distance from the 
observer to the shower axis, and where the retarded dis- 
tance can be rewritten as 



V = ^{-cPst + hy + (1 - n2/32)d2 



(7) 



The retarded time is defined by c{t — tr)/n = i?, where 
R = \J (zP- + n is the index of refraction of air and can 



be expressed as @, Q , 
ct — n^jish 



nV 



(1 - n2/?2) 



(8) 



At a large distance from the shower axis the point-like 
approximation [2], [l^ is valid in which the pancake thick- 
ness is set to zero, fp{h) = S{h). In this limit a simple an- 
alytical expression for the electric field can be derived Q , 



E,it,d) w Jo 



cHH d 



[uMu)] 



(9) 



d^ dtr 

with ctr « —d^/2ct. In general this approximation is valid 
at typical distances above d — 500 m since the typical time 
window for the pulse lies between 10~® and 10~^ s. From 
Eq. © we see that the electric field scales as d~'^ and that 
the integral over time vanishes (i.e. zero response at zero 
frequency) when the shower profile vanishes at the Earth's 
surface. 
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3. Radial dependence of the pulse strength 



Now using 



For tlie sliape of tlie pulse it is very important to under- 
stand the interplay between the different length scales [l^] • 
While at large distances from the shower the important 
length scale is related to a projection of the shower profile 
along the line of sight of the observer, we will show in this 
section that for small distances the pancake thickness pa- 
rameter, L, is important. To describe the electric field for 
small impact parameters we encounter the problem that 
the integral in Eq. (jH)) becomes numerically unstable. We 
therefore first derive a new equation for the pulse which is 
particularly suited at small distances from the core. 

For an observer placed close to the shower axis it is 
more transparent to change variables in Eq. (|6]) and inte- 
grate over the retarded shower time at which the signal is 
emitted since the signal moves with the same velocity as 
the shower. The vector potential can thus be written as, 



A^t,d) 



Jo 



''d{cU) V 



C(t) 

-Jo 
t7{t) 



(10) 



The retarded time U is related to the observing time t and 
the distance from the shower front h by Eq. ([5]) which can 
be rewritten in the limit n = 1 and /3s = 1 as 

= 2(ct + ;.) ' 

The Jacobian in this case is non-zero due to the fact that 
at a fixed observer time the received signal is emitted 
at different shower times, depending on the emission point 
within the pancake. Hence, t~ = — oo and = tr{t,h = 
0) for a fixed observer time t. The Jacobian can now be 
derived from Eq. (jlip . 

dh _ V 

d{ctr) Ctr — h 

The electric field now becomes 



(12) 



E''{t,d) 



-Jo j dUMU)j^ 

t-{t) 

-Jo J duMu)(± + ^^ 



fpih) 



t-(t) 
tt(t) 

- Jo j dtr 



dftiU) fpjh) 



.ft{tr)fp{h) 
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(13) 
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d 
di 

we can write 



_d_ 

dtr 



I3 dfp{h) 
z dh 



(14) 



E^itJ) = -Jo f dtrftitr) 



z dh 



tr{t) 
tt(t) 

Jn / dt 



dftjtr) fpjh) 
dtr Z 



(15) 



where the terms coming from the derivative working on 
the integration limits and the partial integration vanish 
due to the fact that the fp vanishes at the limits, and 
where t, tr and h are related by Eq. ([5|). For very small 
impact parameters the second term on the right-hand side 
in Eq. ((T5|) becomes an integral over the time derivative of 
the shower profile and vanishes at the given limits. This 
happens since the signal travels along with the shower and 
thus h becomes approximately independent of tr- Hence 
the main contribution to the electric field for very small 
impact parameters comes from the first term on the right- 
hand side of equation Eq. ([T5]) . 

When the shower profile is finite at the surface of the 
Earth, z = 0, the terms in Eq. (jlSp suffer from a 1/z di- 
vergence. To regularize this divergence we introduce an 
exponential suppression term in the shower profile func- 
tion. 



Ftitr) = Mtr) - MO)e 



-z/a 



(16) 



where a is the cut-off distance. The divergence in Eq. (jlSp 
is now eliminated since the shower profile function van- 
ishes linearly for z = 0. The electric field for different 
values of a is given in Fig. [T] where the simulations are 
done for a vertical incoming shower with a primary energy 
of Ep = 10^® eV. The calculations show that at a distance 
of d = 50 m even a value as large as a = 1000 m hardly 
affects the pulse. At a distance of only d — 10 m a value 
of a = 100 m gives already a stable result. It has been 
verified that this result hardly depends on the altitude of 
the observer or energy of the cosmic ray, and can also be 
used for more energetic showers where the shower maxi- 
mum is close to the Earth's surface. This is due to the fact 
that the maximum of the observed electric field is emitted 
well before the shower is fully developed. We note here 
that, therefore the observation of radio signals from air 
showers are particularly suited to study the early stages of 
the shower dcvelopement. Interesting to notice is that the 
signal at a distance of only 50 m from the core is hardly 
affected by the shower properties at a height below 1 km. 

Another interesting aspect of Fig. [T] is that the pulse- 
shape docs not vary strongly with distance. At large dis- 
tances where Eq. ^ applies one expects a variation like c?^ 
in pulse length. Instead one finds that the zero-crossing 
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Figure 1: The unfiltcrcd electric field for different values 
of the cut-off parameter a — 10, a — 100, and a = 1000, at 
a distance d ~ 10 m (top) and d = 50 m (bottom) from the 
shower axis. The pancake thickness parameter is taken as 
L = 3.9; see Section [5j 



of the pulse changes from 9 to 1 1 ns while the distance 
changes by a factor 5 from d = 10 to 50 m. This time 
scale corresponds to a length scale of about 3 m which is 
determined by the pancake thickness as supported by test 
calculations. 

The dependence of the calculated pulse height on d is 
shown in Fig. [2l Based on Eq. (|9]) one would expect the 
pulse strength to be proportional to d~^ as given by the 
black line in the figure. Only at very large distances this 
dependence is seen. 

In Fig. [3] the power P above 25 MHz divided by that 
at a distance of d = 50 m is plotted function of d. 
This quantity can easily be compared to actual measure- 
ments since the power is not affected by dispersion of the 
signal contrary to the pulse height and it corresponds to a 
semi-realistic filtering of the pulse. It is clearly seen that 
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Figure 2: Lateral distribution function for the pulse height 
as function of distance on double logarithmic scale for dif- 
ferent values of the cut-off parameter a. 
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Figure 3: Lateral distribution function for the received 
power at frequencies above 25 MHz divided by that at a 
distance of c? = 50 m. 



for small distances a smaller value of a should be used to 
reach a convergent result. The LDF for the power shows 
a distinct exponential structure for distances beyond 5 m 
from the core where the power drops by one order of mag- 
nitude over a distance of close to 70 m. 

For observer distances approaching d = m the pulse 
height tends to diverge. Since, for n = 1, the signal travels 
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along with the shower with the same speed, an observer 
placed at the impact point of the shower will see the equiv- 
alent of a divergent Cherenkov peak. The precise struc- 
ture of the pulse at very short distances will thus strongly 
depend on the index of refraction Q- In addition one ex- 
pects that at close proximity to the core the electric field 
is influenced by the lateral extent of the electrons in the 
shower since the particle density is strongly peaked near 
the shower axis. In a forthcoming publication wc will look 
into this effect in more detail. 



4. Interference 

The main secondary contribution, which is important 
for the present discussion, is due to the charge excess in 
the shower [2, Q which is due to knock out from ambient 
air molecules. The coherent emission from the charge ex- 
cess in cosmic-ray induced showers was first discussed by 
Askaryan (2^, where it is due to the creation and annihi- 
lation of energetic electrons. Even though the Askaryan 
effect is often associated with Cherenkov radiation, it also 
applies to the case in which the index of refraction is close 
to unity 

Macroscopically the charge-excess contribution can be 
described using a net negative charge moving with the 
speed of light toward the surface of the Earth. The con- 
tribution to the vector potential is now given by. 




dt' (17) 



= i^cjMzimm. ,18, 



47reo 



V 



The electric field is now obtained by. 



47reo J R' 



^ dh:^UiU)fih) (20) 



47reo J R' 



^ ' dz del 

OO 



47reo J R 



dh-.-MU)fih) (21) 



where Cx = 0.23 is the approximate fraction of charge 
excess with respect to the total number of electrons and 
positrons in the shower as predicted by the Monte-Carlo 
simulations described in Section [5l 
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Figure 4: The strength of the signal for two different po- 
larization directions is plotted as a function of the orien- 
tation of the observer with respect to the shower. The 
y-polarization has been multiplied by a factor 100 to put 
it on a similar scale. The magnetic field is chosen pointing 
to the North, while (f) = 90° {(j) = 180°) corresponds to the 
observer at the East (North) side of the shower. 

The interference of the gco-magnetic and the charge- 
excess contribution will be important for the observed pulse. 
For this interference the polarization of the two contri- 
butions should be considered. For a vertical shower (i- 
direction) in a horizontal magnetic field (y-direction) the 
electric pulse generated by the induced electric current is 
polarized in the ^-direction independent of the orienta- 
tion of the observer with respect to the shower axis. The 
radiation due to charge excess is polarized radially. De- 
pending on the position of the observer the two contribu- 
tions may interfere destructively (observer on the negative 
i-axis, y = 0), constructively (observer on the positive x- 
axis, ?/ = 0), or be orthogonal (observer anywhere on the 
2/-axis). This effect can be seen clearly in Fig. 2] where the 
dependence of the signal strength in the and y-direction 
is plotted as function of the orientation of the observer. 

In Fig.[S]the LDF is shown separately for geo-magnetic 
and charge-excess radiation showing that while the two 
contributions differ by a factor 5 near the shower core, 
their magnitudes are about equal at a distance of 500 m. 
Depending on the orientation (azimuth) of the observer 
with respect to the shower the two contributions will in- 
terfere constructively or destructively. Since the charge- 
excess and the geo-magnetic contribution depend in a dif- 
ferent way on the shower profile their LDF's differ. This 
gives rise to the interesting observation that for an ob- 
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Figure 5: The behaviour of the LDF as a consequence of 
different types of interference between geo-magnetic radi- 
ation and charge excess radiation. 

server at one side of the shower the charge-excess and 
the geo-magnetic contributions may interfere destructively 
resulting in a LDF that falls off rather steeply close to 
the core and flattens after the point of maximal destruc- 
tive interference where even a local minimum could occur. 
For an observer positioned at the other side of the same 
shower the two contributions add constructively resulting 
in a relatively smooth LDF. In the perpendicular direction 
the two contributions will not interfere as their polariza- 
tions are orthogonal. The distance where a local minimum 
could occur and the LDF starts to flatten will depend on 
the relative magnitude of the charge-excess and the geo- 
magnetic contribution where the latter depends strongly 
on the angle of the shower with the magnetic field. It fol- 
lows that one should be extremely careful in determining 
the LDF [sOl since it is expected to depend on the orien- 
tation of the observer to the shower. 

5. Hybrid approach 

In the previous sections, all simulations were done us- 
ing the basic MGMR model. The particle distributions 
in these simulations were parametrized using the analyti- 
cal formula given in Section [51 In the current section we 
will discuss the hybrid approach to the MGMR model. 
Within this approach the important macroscopic proper- 
tics describing the electric field, such as the pancake thick- 
ness and the shower profile, are obtained from Monte-Carlo 
simulations. 



These simulations are done using the cascade mode of 
the CONEX shower simulation program [sil [s^l- This 
code has been modified to include the deviation of charged 
particles in the Earth's magnetic field, which is discussed 
in detail in Appendix B of Ref. Q . Furthermore an anal- 
ysis tool has been written, giving a histogrammed output 
of the full three-dimensional and timing information of the 
currents and particle distributions within a user-defined 
region of the shower front which in the simplest case be- 
comes an average of all particles over the complete shower 
front. Several examples for such distributions, smoothed 
by fitting procedures, can be found in Ref. 0], more de- 
tails will be discussed in a future publication. The com- 
plete shower-simulation package including analysis tools is 
called CONEX-MC-GEO. 
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Figure 6: Total number of electrons and positrons as func- 
tion of shower depth for 40 proton (blue curves) and 40 
iron (red curves) initiated showers at energies E = 10^^ eV 
(top) and E = lO^^eV (bottom). 

As a first example the histogrammed output obtained 
by CONEX-MC-GEO, which will be used in Section |6] to 
study the differences of the LDF for a proton- and an 
iron-primary particle, will be discussed. The simulations 
arc done for 40 proton- and 40 iron-induced showers to 
extract the longitudinal shower profiles and the pancake 
thicknesses at energies of 10^^ and 10^® eV. 

The obtained shower profiles are shown in Fig. [51 As 
is expected the shower-to-shower fluctuations are rather 
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large for proton- and small for iron-induced showers. It 
should be noted that the intrinsic fluctuations for a few 
showers become large. This is partially due to thinning 
effects. The Xmax is larger for proton than for iron. The 
shower profile drops to zero at around 870 g cm^^ where 
the shower hits the surface of the Earth according to the 
conditions at the site of the Pierre Auger Observatory. 

In the previous section the importance of the charge- 
excess contribution and its impact on the LDF is shown. 
To make this more quantitative the fraction of excess elec- 
trons is plotted in Fig. [7] for Monte-Carlo simulations of 
40 proton- and 40 iron-induced showers at an energy of 
10^^ eV. An analysis for 10^® eV showers gives similar re- 
sults. The charge excess is slightly increasing with increas- 
ing shower depth. Near the shower maximum the charge 
excess is close to 23% of the total number of electrons and 
positrons which is the value we have assumed in all our cal- 
culations described in Section The charge-excess con- 
tribution scales linearly with the charge-excess fraction. 
The shower-to-shower fluctuations thus have a relatively 
small effect. The other important parameter for the radio- 
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Figure 7: The fraction of charge-excess electrons with re- 
spect to the total number of electrons and positrons as a 
function of shower depth for 40 proton and 40 iron induced 
showers at an energy of 10^^ cV. 

emission calculations is the pancake thickness parameter. 
The relevant parameter for this is the mean distance of the 
electrons behind the shower front, (h), the mean pancake 
thickness parameter. In Fig. [5] this quantity is plotted for 
the set of showers at 10^^ eV. Again the fluctuations are 
considerably larger for protons than for iron. In addition 
one can see that the mean pancake thickness parameter is 
almost independent of shower height and of shower type. 
Calculations for showers for 10^^ eV give similar results. 

In the present calculations L = (h) has been fixed 
for the full shower development. From Fig. |8] it can be 
seen that L = 4.3 m is a reasonable average value for iron- 
induced air showers, and L ~ 3.9 for proton-induced air 
showers. Where the width of the function at 1/e of the 
maximum can be shown to be approximately 1.5L. The 
fluctuations in the average value of L are small with re- 



spect to the fluctuations in the shower profiles and not 
considered in further analysis. 



7.0 




-25 -20 -15 -10 -5 

t(/iS) 



Figure 8; Mean distances of electrons behind the shower 
front for iron (red curves), and proton (blue curves) in- 
duced showers, as function of the negative shower time. 



6. Composition 

As a first application of the hybrid approach we study 
the LDF for proton- and iron-induced showers using MGMR 
in combination with CONEX-MC-GEO. The smoothed 
shower profiles extracted from each Monte-Carlo simula- 
tion are used in the calculation of the radio signal. For the 
calculation at short distances we use a suppression coefB- 
cient of a = 100, since in Section |3] it is shown that this 
value gives reliable results at the distances considered here. 
The maximum of the field strength as function of distance 
is shown in Fig. [HI This is done for an observer placed 
on the X-axis where the charge-excess and geo-magnetic 
fields interfere constructively. This figure shows that the 
slope for iron-induced showers is less steep than for pro- 
tons. One also sees that at a distance of d = 100 m the 
calculated signal strength for iron and proton showers is 
almost the same. These features were also observed in 
REAS2 calculations 0. 

To express the differences in the LDF more clearly we 
introduce a variable that is relatively easy to extract from 
data, 

- ^(50, / > 25)/P(300, / > 25) (22) 

where P{d, f > 25) is the power in the pulse for frequen- 
cies larger than 25 MHz at a distance dm from the shower 
core. In Fig. [10] the results of the calculations are dis- 
played in a histogram. The averages of the values for iron- 
and proton- induced showers differ strongly, however due to 
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Figure 9: The LDFs are shown for proton (blue diamonds) 
and iron (red crosses) induced showers at an energy of 
lO^'^eV (top), and lO^^eV (bottom), displayed in Fig. [SI 
The LDFs are given for an observer position leading to 
constructive interference between the charge excess con- 
tribution and the contribution due to gco-magnetic radia- 
tion. 



in Fig.ini As noted earlier the LDF, and thus also R^q/^qq, 
depends on the orientation of the observer with respect 
to the shower and comparisons should thus be made at 
the same azimuth with respect to the shower core. It has 
been checked that the qualitative differences between iron 
and protons remain as shown in Fig. [10] also for other az- 
imuth positions of the observer. The reason for this is 
that the charge-excess fraction is similar for proton- and 
iron-induced showers. 
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Figure 10: .R50/300J see text, for the 40 simulated proton 
and iron showers at an energy of i? = 10^^ eV (top) and 
E = IQi^eV (bottom). 



shower-to-shower fluctuations some of the proton showers 
give rather similar values as for iron. The proton showers 
that give a similar value for ^^50/300 a-re showers which have 

profiles (see Fig. ^ very similar to those of iron. At higher 7_ Summary and Conclusions 

energies the average value of ^59/390 differs less between 



proton- and iron-induced showers. This is a direct con- 
sequence of the fact that Xmax-, the important parameter 
responsible for the observed effect, differs less for the dif- 
ferent primary particles at higher energies, as can be seen 



We have derived a new expression for the electric field 
emitted by an EAS at close proximity to the shower axis 
in the MGMR model. It differs from the expression used 
for the pulse strength at large distances by a change of 
integration variables. This change of variables expresses 
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that at large distances the pulse shape is basically deter- 
mined by the shower profile where the pancake thickness is 
integrated over. At short distances the picture is reversed 
and the pulse shape is determined by the pancake function 
and the shower profile can be integrated. 

Due to the fact that the electron density is finite when 
the shower reaches Earth the pulse height becomes di- 
vergent at short distances from the core. This we have 
resolved by -artificially- suppressing the shower at small 
heights above the surface of the Earth. We show that this 
allows us to get accurate numerically-stable results at dis- 
tances down to 30 m. At smaller distances other effects 
such as the lateral spread of the electrons in the shower, 
as well as shower-to-shower fluctuations will be important 
which are not included in the present schematic approach. 
In a future publication these effects will be considered in 
detail. 

Since at short distances, the pulse shape is determined 
by the pancake function while at large distances it is gov- 
erned by the shower profile this offers the possibility to 
use the LDP of the radio pulse to distinguish iron- and 
proton-induced showers. Before addressing this point we 
have added the radiation due to charge excess as the most 
important secondary process. It is shown that the determi- 
nation of the polarization as a function of the azimuth an- 
gle can be used to study the different emission mechanisms 
of radio emission from an EAS. We show that the LDF is a 
very powerful tool to distinguish iron- and proton-induced 
showers although the effects of shower-to-shower fluctua- 
tions for protons can be large, where at the same time the 
power at a distance of 50 m from the shower core can be 
used as a measure for the energy. 
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